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We investigate the quantum phases of monodispersed bosonic gases confined to a triangular lattice 
and interacting via a class of soft-shoulder potentials. The latter correspond to soft-core potentials 
with an additional hard-core onsite interaction. Using exact quantum Monte Carlo simulations, we 
show that the low temperature phases for weak and strong interactions following a temperature 
quench are a homogeneous superfluid and a glass, respectively. The latter is an insulating phase 
characterized by inhomogeneity in the density distribution and structural disorder. Remarkably, 
we find that for intermediate interaction strengths a superglass occurs in an extended region of the 
phase diagram, where glassy behavior coexists with a sizable finite superfluid fraction. This glass 
phase is obtained in the absence of geometrical frustration or external disorder and is a result of the 
competition of quantum fluctuations and cluster formation in the corresponding classical ground 
state. For high enough temperature, the glass and superglass turn into a floating stripe solid and 
a supersolid, respectively. Given the simplicity and generality of the model, these phases should be 
directly relevant for state-of-the-art experiments with Rydberg-dressed atoms in optical lattices. 
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It is well established that bosonic and fermionic sys¬ 
tems subjected to a disordered external potential feature 
localization phenomena [1, 2]. The interplay between dis¬ 
order, interactions and many-body quantum effects such 
as superfluidity is now a subject of intense research [3- 
8], as, e.g., bosons in random environments occur in a 
variety of experimentally relevant systems ranging from 
cold atoms [9-13], to superconductors [14] and quantum 
liquids [15, 16]. Usually, the combination of disorder and 
repulsive interactions inhibits the emergence of superflu¬ 
idity and Bose-Einstein condensation (BEG) and leads to 
an insulating gapless phase, known as Bose glass [17-19]. 

Remarkably, results of quenched Monte Carlo simula¬ 
tions in the context of “^He have shown that superfluidity 
and BEG may coexist with structural disorder and in¬ 
homogeneity (i.e., glassy physics) in the absence of any 
random external potentials [20]. The resulting out-of- 
equilibrium state was termed superglass (SG), as a dis¬ 
ordered analog of the supersolid (SS) phase [21]. While 
experiments have so far remained inconclusive [22, 23], 
this proposal has spurred considerable theoretical activ¬ 
ity to derive possible microscopic models of a SG [24-30]. 
Exact numerical results for bosons on lattices have shown 
that a thermodynamic SG phase can indeed emerge as a 
result of a competition of quantum fluctuations and ex¬ 
ternally induced frustration. For attractive interactions 
the latter can be induced via a random chemical potential 
[27], while for repulsive ones a SG can occur in theoret¬ 
ical models where either a self-disordered environment 
is induced by geometrical frustration (e.g., on random 
graphs) [28] or where disorder is a consequence of prop¬ 
erly chosen random inter-particle interactions [29, 30]. 
In this context, main open questions are whether it is 
possible to obtain a SG in any theoretical models where 
frustration is not artificially built in the Hamiltonian, 


and if this new phase of matter may be experimentally 
observable in any physical system. 

Here, we show that the SG phase can exist for a large 
class of bosonic lattice Hamiltonians. The latter are of 
the extended Hubbard-type, featuring a soft-shoulder in¬ 
teraction potential. Surprisingly, glassy behavior is ob¬ 
tained in the absence of any externally imposed frustra¬ 
tion e.g., in the lattice geometry, or in the interactions. 
Rather, frustration is here induced by cluster formation 
for large particle density, similar to the conditions of SS 
formation in soft-core models [31-33]. As an example, we 
consider a simple triangular lattice with isotropic two- 
body interactions. We analyze the phases and, follow¬ 
ing a quench in the temperature T or in the interaction 
strength, demonstrate the existence of both a classical 
glass (G) and a SG at low enough T. The latter are 
the out-of-equilibrium counterparts of a floating stripe 
solid (S) and a SS, respectively. These glass and super¬ 
glass phases should be observable in experiments with 
Rydberg-dressed alkali atoms loaded into optical lattices. 

The relevant Hamiltonian for hard-core bosons con¬ 
fined to a 2D triangular lattice reads 

H = -tY^{h^ihj+b]b^ FV ^ n,nj. ( 1 ) 

{i,j} i<r, rij<r„ 

Here, bi (6j) are hard-core bosonic annihilation (creation) 
operators at site i, rii = blbi, rij is the distance between 
sites i and j, and t is the tunneling rate on a lattice of 
spacing a. In the following, t and a are used as units 
of energy and length, respectively. In classical physics, 
the soft-shoulder potential of Eq. (1) is of interest for 
soft-matter models of, e.g., colloids [34-36]. In quantum 
physics, this potential can be engineered in clouds of cold 
Rydberg atoms, where both the strength V and the range 
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Tc of the interaction can be tuned by weakly-admixing the 
Rydberg level to the ground state [32, 33, 37-42, 44] (see 
SupMat). The additional onsite hard-core constraint can 
be enforced using, e.g., Feshbach resonances. 

The quantum phases of Eq. (1) with Tc = 1 (i.e., 
nearest-neighbor interactions) are well known [45-49]: 
for densities p < 1/3 (p > 2/3), p = 1/3 (p = 2/3) and 
p > 1/3 (p < 2/3) the low-energy phase is a superfluid 
(SF), a gapped lattice S, or a gapless SS, respectively. 
The latter is an exotic state of matter where density cor¬ 
relations (here with \/3 x -\/3 ordering) coexist with a 
finite superfluid fraction ps, which is a result of doping 
the solid with interstitials (vacancies). The SS phase is 
generally robust against perturbations to the Hamilto¬ 
nian (1), and may be observed experimentally, e.g., with 
cold quantum gases trapped in optical lattices and inter¬ 
acting via dipolar interactions [50-52]. 

In this work, we are interested in Eq. (1) with Tc > 1. 
For Tc > 1 the interactions belong to a large class of po¬ 
tentials that support the formation of self-assembled clus¬ 
ters of particles for sufficiently large densities Vc^/p > 1 
[35, 36]. Such a phenomenon is essentially independent of 
the details of the interactions, as long as the latter display 
a negative Fourier component [34]. In the classical regime 
(i.e., t = 0) cluster formation has been shown to lead to 
frustration, which is manifested in an exponential growth 
of the ground state degeneracy as a function of the sys¬ 
tem size [53]. In the quantum regime (i.e., t > 0) this 
leads to several novel exotic phenomena at equilibrium: 
anomalous Luttinger-Liquid behavior [53] and emergent 
supersymmetry in ID lattice geometry [54] as well as free- 
space supersolidity in 2D [32, 33, 55]. The latter occurs, 
for appropriate values of interaction strength, at any den¬ 
sity fulfilling the clusterization condition rc^/p > 1 [55]. 
In the following we consider, as a way of example, the 
simplest cluster forming potential with Tc = 2 and in¬ 
commensurate particle densities consistent with such a 
condition. Our main focus is the demonstration of a G 
and SG emerging when a crystal and a SS are driven out 
of equilibrium via a temperature quench, respectively. 
Glassy phases for different p, and quench protocols 
are discussed in the Supplementary Material (SM). 

We study the Hamiltonian in Eq. (1) by means of Path 
Integral Quantum Monte Carlo simulations based on the 
worm-algorithm [56]. This technique is numerically ex¬ 
act for bosonic systems and allows for accurate estimates 
of the superfluid fraction ps = ((IF| -I- IF^)/(6/3p)) and 
the static structure factor S{k)/N = - exp[— ik • (r^ — 

Vj)]{ninj)/N'^ {N is the lattice size). The latter mea¬ 
sure superfluidity and diagonal crystalline orders, respec¬ 
tively. Here, 13 = 1/(A:bT), with ks Boltzmann con¬ 
stant (in the following set to 1), H/ is the winding num¬ 
ber in the *-th direction, k is a lattice wave-vector, and 
(• • •) stands for statistical average. In addition, we com¬ 
pute the renormalized Edwards-Anderson order parame¬ 
ter qEA = (IeaIoas^-i which, in the absence of crystalline 



FIG. 1: (color online), (a): Superliuid fraction pa, and renor¬ 
malized Edwards-Anderson parameter qEA as a function of 
V/t, for T/t = 1/12. (b): ps and qea as a function of T/t, 
for V/t = 5.4. In both panels the density is p = 13/36 and 
the lattice size is A = 900. Solid lines are guides to the eye. 
Inset: maximum value of the structure factor Suiax/N as a 
function of l/y/N for p = 13/36, V/t = 5.4 and T/t = 1/12; 
the dashed line is a linear fit for the three largest system sizes. 


order, is the well-accepted observable to identify glassy 
behavior on a lattice [28, 57]. Here, qEA = P)"^ 

and = N p{l — p) is its maximum value obtained for 
a classical situation with no particle delocalization. We 
perform large-scale simulations with up to A = 2304 lat¬ 
tice sites and temperatures as low as T/t = 1/12. For 
each N and T, numerical values for the observables above 
are obtained by averaging over a minimum of 32 and a 
maximum of 100 different realizations of the quench. 

Figure 1 [panel (a)] shows example results for the 
superfluid fraction ps and the renormalized Edwards- 
Anderson parameter qEA as a function of the interaction 
strength V/t for N = 900 and T/t = 1/12. Within the 
interesting range of interaction (5.0 < V/t < 6.0), ps 
is found to decrease monotonically with increasing V/t 
from approximately 0.25 to about 0.05. In the same pa¬ 
rameter range, qEA increases up to values of the order 
of ^ 0.2. We note that in this regime the system does 
not feature crystalline order, i.e., the computed structure 
factor S'(k)/A vanishes for any non trivial wave vector 
k 0 in the thermodynamic limit, as proven in the Inset, 
where we show the scaling with of 5'max/A, the 

average of the largest peaks of the structure factor over 
several quench realizations. In the same limit the super¬ 
fluid fraction stays finite. These data demonstrate one of 









3 


5.00 



FIG. 2: (color online) Phase diagram of model Eq. (1) with 
Tc = 2 as a function of temperature T/t and interaction 
strength V/t, for particle density p = 13/36. Equilibrium 
phases: normal liquid (L), superfluid (SF), stripe-crystal (S), 
and supersolid (SS). A temperature quench to final values of 
T/t below the dashed line leads to a glass (G). The existence 
of a superglass (SG) is demonstrated below the dotted line 
using the same quench protocol. 


the main results of this work, namely the existence, in an 
extended region of parameters, of a SG, corresponding to 
an inhomogeneous non-crystalline superfluid. The depen¬ 
dence on T/t of both qea and ps is shown in Fig. 1 [panel 
(b)]. For the specific value of the interaction strength in 
the panel superglassiness is realised below T/t ~ 0.2. 

For weak interactions, the SG phase quantum melts 
into a regular homogeneous superfluid (SF) with Ps > 0 
and Pea — 5'(k) = 0. For the parameters of Fig. 1(a) this 
is obtained by decreasing the interaction strength below 
V/t ~ 4.8. On the other hand, sufficiently large interac¬ 
tion strengths are found to inhibit superfluidity and turn 
the SG into an insulating G. The latter is characterized 
by a finite value oiqEA and ps ~ S'(k) = 0 [i.e., V/t >6.2 
in the figure]. Within this glass phase quantum effects 
are largely suppressed. While classical glasses are well 
known to appear in disordered spin models, as well as in 
certain polydispersed systems of particles [57], here we 
demonstrate that (classical) glassy physics may emerge 
in a simple and rather general model of immediate ex¬ 
perimental interest for bosons on a regular lattice. 

The computed phase diagram of Eq. (1) is shown in 
Fig. 2 for a choice of particle density p = 13/36 as a func¬ 
tion of T /t and V/t. At high temperatures we find a nor¬ 
mal liquid (L) phase independently of the values of F/t, 
as expected. For sufficiently small interaction strength 
V/t < 4.8, this normal phase turns into a homogeneous 
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FIG. 3: (color online) (a): Normalized Edwards-Anderson 
order parameter qea for model Eq. (1) as a function of 
temperature T. The value of the interaction strength is 
V/t = 10. Data are for a lattice with size N = 900 and 
density p = 13/36. The solid line is a guide to the eye. In¬ 
set: Size dependence of Smax/N for V/t=10. Values of the 
temperature are T/t = 1/12 and T/t = 3/2 for the full and 
empty symbols, respectively. Dashed lines are fits to our nu¬ 
merical data, (b-c): Maximum value of the structure factor 
SinsL/N obtained in a given realization of a quench, plotted 
as a function of the number of different quench realizations. 
The latter only differ in the (random) initial condition and 
in the thermalization seed. The corresponding phases in the 
thermodynamic limit are indicated. Fluctuations in the val¬ 
ues of S^:x/N indicate glassy behaviour. The parameters are: 
N = 1764, T/t = 1/2.5, V/t = 5 (SS), N = 2304, T/t = 1/9, 
V/t = 5.4 (SG), N = 1296, T/t = 1/0.7, V/t = 10 (S), 
N = 1296, T/t = 1/12, V/t = 10 (G). 


superfluid by decreasing T /t, via a phase transition which 
is consistent with the Berezinskii-Kosterlitz-Thouless sce¬ 
nario. On the other hand, for large enough V/t and fol¬ 
lowing a quench to low T the system displays a marked 
insulating glassy behavior with pea 0, ^(k) = ps = 0 
[Fig. 3(a) and full symbols in the Inset]. The interplay 
between glassy physics and superfluidity is mostly evi¬ 
dent for values of T/t below the dotted lines in Fig. 2, 
resulting in the SG scenario discussed above. 

Interestingly, we find that fluctuations can restore crys¬ 
talline order for sufficiently large T. This is shown for in¬ 
termediate temperatures in Fig. 2, where a S (SS) phase 
intervenes between the low-temperature G (SG) and the 
high-temperature normal L. Here, the crystal is a floating 
stripe solid, with finite diagonal long range order in the 
thermodynamic limit. Examples for the finite size scaling 
of the maximum value of the structure factor Smax/N in 
the S and G phases are shown in the inset of Eig. 3(a) 
(empty and full symbols, respectively). While in the S 
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phase Smax/N is essentially independent of the system 
sizes investigated in this work, in the G phase S^nax/N 
vanishes in the thermodynamic limit. In both phases 
Ps - 0. 

The difference between the glassy and crystalline 
phases is shown in Fig. 3(c), where we plot the maxi- 

( D\ 

mum value of the structure factor S'max/iV for each indi¬ 
vidual realization of a temperature quench at a given V/t 
and final T/t. In the crystalline phases, S^^x/N is es¬ 
sentially identical in all realizations and S^ax/N remains 
finite in the thermodynamic limit. However, within the 
glassy phases, S^x/N can fluctuate widely and in aver¬ 
age decreases to zero with the system size. As shown in 
Fig. 3(b) the dependence of Smax/N on the realization 
for the SG and the SS phase is similar to that for the G 
and the S ones, respectively. 

Further insight into the phases of Hamiltonian (1) is 
given by the analysis of the averaged site-density maps 
in Fig. 4. Specifically, we show results for a portion of 
the system and for a choice of T/t and V/t such that the 
system is a SF [panel (a)], a SG [panel (b)], and a G [pan¬ 
els (c-e)]. For comparison, panel (f) shows a cluster-type 
crystalline phase [i.e., 5'(k) yt 0] stabilizable at a density 
p = 1/3, for V/t = 10 and T/t = 1 [58]. In the homo¬ 
geneous SF the average occupation number at each site 
equals the density p of the system, as expected. The re¬ 
sulting value of qEA is thus negligible. Gonversely, when 
V/t is large [panel (c)] the spatial density is highly inho¬ 
mogeneous: particles form self-assembled clusters char¬ 
acterized by different numbers of constituents and spatial 
orientations, as well as by varying inter-cluster distances. 
These features lead, in the thermodynamic limit, to the 
absence of diagonal long range order, similar to an (emer¬ 
gent) polydispersity. Noticeably, the occupation num¬ 
ber of lattice sites between clusters is here substantially 
suppressed, signaling particle localization. The resulting 
glass phase is insulating, similarly to, e.g, a regular Bose 
glass obtained by externally induced disorder [17]. 

Figure 4(b) shows that cluster formation (i.e., inho¬ 
mogeneity) persists even at intermediate values of V/t, 
leading to a non zero value of qsA in the absence of crys¬ 
talline order. The occupation of inter-cluster lattice sites 
is here enhanced with respect to panel (c). Such an en¬ 
hancement is due to the presence of quantum fluctuations 
and exchanges of identical particles, responsible for the 
finite value of ps and thus of superglassiness. 

While here we have mainly focused on the minimal 
model Hamiltonian with = 2 and a given density 
p — 13/36, we have verified that glassy phases also oc¬ 
cur for larger Vc and densities p satisfying the clusteriza¬ 
tion condition, as well as when the details of the cluster 
forming potential are changed (see Fig. 4 and SM). In 
particular, glassy phases are obtained for soft-core van 
der Waals type interactions relevant to experiments with 
Rydberg dressed atoms [42, 43] [see Fig. 4 and SM]. In 



FIG. 4: (color online) Averaged site density for a portion 
of the system. Black circles depict the lattice sites. Den¬ 
sity values are proportional to the size of red circles. Panel 
(a) shows an homogeneous superfluid phase (T/t =1/9 and 
V/t = 4); panels (b) and (c), refer to a SG (T/t = 1/12, 
V/t = 5.4), and to a normal G (T/t = 1/12, V/t = 10), re¬ 
spectively. In panels (a-c) the density is p = 13/36. Panels (d) 
and (e) show the glassy density map obtained for p = 0.401 
(V/t = 10, T/t = 1/12), and for the same density of panels 
(a-c), nsing the van der Waals purely repulsive soft poten¬ 
tial with Vryd/t = 30, T/t = 1/3 (see details in SM), re¬ 
spectively. Panel (f): a crystalline structure stabilizable at 
p = 1/3, V/t = 10 and T/t = 1. 

view of possible experiments, where V/t can be easily 
varied, we have verified that [see SM] (i) quenches inV/t 
for a fixed low T from a superfluid lead to the same SG 
and G phases described above; in addition, (ii) the SG 
and G phases are robust against density variations at the 
percent level. In particular, the G phase is found at all 
densities explored, up to p ~ 0.4, for sufficiently large 
V/t. We emphasize that while the finite lifetime of a 
Rydberg-dressed gas can be a limitation to the observ¬ 
ability of phases in thermal equilibrium, glassy phases 
are realized out of equilibrium and their relevant prop¬ 
erties are essentially those emerging in the first stages of 
the experiment. 

The glass and superglass phases discussed here are 
the low-T quenched counterparts of the equilibrium solid 
and supersolid found for intermediate T. These latter 
equilibrium phases are analogous to those found in, e.g.. 
Ref. [55] in continuous space. While recent results point 
to the existence of out-of-equilibrium glassy-type phe¬ 
nomena in the classical regime for cluster-forming sys¬ 
tems [59] (see also Ref. [60] for a discussion on the dy- 
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namics of the equilibrium phases), the existence of glassy 
dynamics in the quantum regime and in continuous space 
remains an open question. 

In conclusion, we have demonstrated that glassy 
phases can be realized for a broad class of simple 
bosonic frustration-free Hamiltonians of the extended 
Hubbard-type. For intermediate interaction strength 
the interplay between quantum fluctuations, statistics 
and glassy physics gives rise to an exotic SG scenario, 
where glassiness coexists with superfluidity, in contrast 
to a conventional Bose glass. In our model, frustration 
arises from the self-assembling of clusters, which is 
a direct consequence of the (isotropic) inter-particle 
interaction potential at high enough density. The physics 
described in this work should be directly relevant for 
experiments with ultracold Rydberg-dressed atoms 
confined to optical lattices [32, 41, 61]. We hope that 
our work will provide new insights for unveiling other 
general mechanisms to glassy physics, and in general to 
frustration-induced phenomena both in the classical and 
quantum regime. Interesting extensions might include 
the search for exotic phenomena beyond the SG, such as 
frustration-induced Bose metals [62, 63] and emergent 
gauge fields [64]. 
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We demonstrate the occurrence of glassy physics in the model Hamiltonian introduced in Eq. (1) of the main text for different 
choices of interaction radius Cc > 1 and particle density p fulfilling the clusterization condition r^y/]} > 1. We find no evidence 
of a glass for = 1, i.e., when cluster formation does not occur. These results indicate that glassiness is favored by the 
formation of self-assembled clusters. Low-temperature glassy scenarios are found following either a quench in temperature T 
or in the interaction strength Vjt. Glassy phases are robust against variations of the details of the cluster-forming interaction 
potential. Parameters for a possible realization of glassy dynamics with Rydberg-dressed gases trapped in optical lattices are 
discussed. 


Dependence of glassy behavior on the interaction 
radius rc, density p and cluster formation 

In this section we investigate the dependence of the 
glass behaviour found in the main text on the radius 
of the interaction potential of Hamiltonian (1) and on 
the particle density p. In particular, we are interested in 
exploring whether the appearance of glassiness is related 
to the formation of clusters in the ground state. 

By performing extensive quantum Monte-Carlo calcu¬ 
lations, we demonstrate that glassy phases also occur for 
values of Tc and p larger than those used in the main text, 
when the clusterization condition Tcyfp > 1 is met. In 
contrast, for the largest interaction strength and smallest 


temperature investigated in the main text glass behavior 
is not found by choosing Tc = 1, at which value self- 
assembled clusters are absent. These results underline 
the importance of cluster formation for the emergence of 
glassiness in model (I). 

Glass behavior for Vcy/p >1: increasing r^ 

Figure Sl(a) shows results for the maximum value of 

/ D\ 

the structure factor SmA/N for different realizations of 
a low-temperature quench for interaction radius = 3 
(larger than that in the main text) with particle den¬ 
sity p = 13/36 (same of that in the main text). The 
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interaction strength is chosen Vjt = 10, the temperature 
T/t = 1/12, and the system size N = 1296. 

We find that, despite the relatively small value of 

{ T>\ 

N, Sinsix/N clearly fluctuates between different realiza¬ 
tions [Fig. SI (a)]. In addition, panel (b) shows that the 
structure factor averaged over the individual realizations 
Smax/1^ vanishes in the thermodinamic limit. In the 
same limit, the superfluid fraction is ps — 0, while the 
(realization averaged) Edwards-Anderson parameter (jea 
stays finite. These results demonstrate the existence of a 
glass phase with Tc = 3. Here we have utilized the same 
parameters of the glass phase for = 2 in the main text, 
showing that in this case an increase of Vc does not alter 
the physical scenario. 


Glass behavior for Vc^/p >1: increasing p 

Figure S2 shows values of the Ed wards-Anderson pa¬ 
rameter and S^-x/N for different realizations of a 
temperature quench with target T/t = 1/12, V/t = 5.4, 
Tc = 2 (as in the main text) and p = 0.3650. The 
latter corresponds to a variation at the percent level 
of the mostly used density value in the main text (i.e., 
p = 13/36). The system size is A^ = 2304. 

Also in this parameter regime the maximum value of 
the structure factor [panel (b)] depends on the quench 
realization, showing the same glassy behaviour we 
observed at p = 13/36 [see main text, Figure 3, panels 
(b-c)]. The fluctuations of are much less pronounced 
[panel (a)] and the values remains large. In addition, 
we find that here structural disorder and inhomogeneity 
coexist with finite superfluid fraction demonstrating the 
existence of a superglass. 

At sufficiently large V/t, the existence of glassy phases 
for the cluster-forming potential in Eq. (1) of the main 
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FIG. SI: (color online), (a): Maximum value of the structure 
factor S^/x/N for different realizations. The system size is 
N = 1296. (b): size dependence of the average structure 

factor Smax/lV. In both panels p = 13/36, T/t = 1/12 and 
V/t = 10. The hamiltonian of the system is that in Eq. (1) 
of the main text with r^ = 3. 
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FIG. S2: (color online). Edwards-Anderson parameter 
(a) and maximum value of the structure factor S^x/N (b) 
for different realizations. The system size is A = 2304, p = 
0.3650, T/t = 1/12 and V/t = 5.4. For these parameters the 
system is a superglass in the thermodynamic limit (see text). 


text has also been verified for values of the particle 
density fulhlling the clusterization condition Vcy/p > 1 
as high as p 0.4. For this value of p, an example of 
structural disorder and cluster formation is given in Fig. 
4(d) of the main text. 

Given the results above, as Hamiltonian (1) is particle- 
hole symmetric, one may expect that glass behavior could 
be observable for essentially all densities with Vc-s/p > 1, 
Tc = 2 and an appropriate choice of V/t. The inves¬ 
tigation of this point and in particular of the scenarios 
emerging for very large values of ~ \/]V/2 (where the 
interaction range becomes of the order of the whole sys¬ 
tem size) will be the subject of a future work. 


Absence of glassy phases for r^ = 1 (no cluster formation) 


Finally, to underline how cluster formation plays a cru¬ 
cial role for the observation of the glassy scenarios inves- 
tigated here, we show in Fig. S3 estimates of and 
Smax/N for different realizations of a temperature quench 
with target T/t = 1/12, V/t = 10, N = 2304 and = 1. 
With this choice of parameters the interaction potential 
does not support clusterization and the equilibrium phase 
is a supersolid [SI]. 

As shown in figure, the values of and S^x/N are 
essentially identical in all realizations, i.e., the quench is 
ineffective and the equilibrium physics is restored. Con¬ 
versely, at the same temperature, for Tc = 2 and particle 
density satisfying the clusterization condition, we find 
that the quenched counterpart of a supersolid is a su¬ 
perglass (see Fig. 2 main text). This indicates that clus¬ 
ter formation enhances frustration in the system favoring 
glassy behavior. 
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FIG. S3: (color online). Same observables as in Fig. S2 for 
T/t = 1/12, Vjt = 10, N = 2304 and — 1. In this case 
cluster formation does not occur. As opposed to the = 2 
case where cluster formation induces disorder and result in 
an out-of-equilibrium superglass, here a temperature quench 
does not alter the equilibrium supersolid phase. 


Glasses with soft-shoulder repulsive van-der-Waals 
interactions 


In general the phases studied in this work should be 
relevant for the whole large class of interaction potentials 
leading to cluster formation in the classical ground state 
[S2]. Among these potentials, an interesting example is 
Kyd(?') = 14yd/[l + {i'/t’cT], which is realizable with cold 
Rydberg atoms [S3-S7], where the ground state of each 
atom is off-resonantly coupled to an excited Rydberg 
state using a laser with (effective) Rabi frequency 
and red detuning A, with A [S8-S15]. In the 

appropriate parameter regime, I/.yd = r2^/(167rA^) 
and Tc = (C'6/2A)^/® is usually of the order of a few 
/im, with Cq the coefficient of van-der-Waals type 
interactions for the atoms. We have checked numerically 
that glassy phases similar to those described above 
can be obtained for I4yd sufficiently large. As an 
example. Fig. 4(e) (main text) shows results for the 
density within the glass phase with Id-yd/t = 30. Below 
we give an example of possible parameters to achieve 
the regime of strong interactions in an optical lattice. 
While many more schemes for Rydberg dressing are 
currently being developed [SIS, S14, SI6], we notice 
that hrst examples of Rydberg-dressing have been very 
recently realized in three different experiments [SI7-S19]. 

Here we consider a gas of ®^Rb atoms, where the 
groundstate Ig) of each atom is coupled off-resonantly to 
an excited Rydberg state |r) using a (two-photon) laser 
with (effective) Rabi frequency H and detuning A from 
the atomic transition. We choose red detuning in order 
to realise the wanted soft-shoulder potential [SIO, S20] 
and |A| ^ H, in order to satisfy the weak dressing 
condition. As an example, we consider |r) = 14351/2), 
r2/(27r) = 2.0MHz and A/(27r) = 20 MHz, where we 
have set Planck’s constant h= 1. These parameters im¬ 


ply Vryd — 450 Hz. The resulting cut-off radius Tc > 2a, 
with a ~ 0.5^m the characteristic length for the lattice 
spacing, allows one to easily access the cluster regime. 
The hopping amplitude t depends exponentially on the 
depth of the optical lattice [S21], which allows one to 
achieve the regime of strong interactions V/t. For exam¬ 
ple, for an in-plane lattice depth Vq/Eji > 13, Eji = 2ttx 
50 kHz, tl{27r) < 40 Hz, and thus V/t > 10. This 
value can be further increased by increasing Vq- A depth 
Vj_ of the optical lattice in the transverse direction with 
Vj_ ^ Vq ensures that the dynamics is two-dimensional. 
Here, i?fl/(27r) = 7r^/(2ma^) is the recoil energy of atoms 
with mass m. Finally, the small effective decay rate 
7 eff/( 27 r) = (H/A)^7j. of the dressed ground-state re¬ 
sults in a 7 eff usually of the order of tens of Hz, with 
7 r the bare decay rate of the Rydberg state. While this 
small rate can be still detrimental for reaching equilib¬ 
rium phases, as noted in the main text the glassy phases 
of interest here are obtained out of equilibrium. The en¬ 
suing dynamics in the glassy phase after the quench is 
essentially frozen and could be observed immediately af¬ 
ter the quench. We further note that, as shown above, 
variations of the particle number by a few percents does 
not hinder the realization of glassy phases.Provided that 
residual spontaneous emission from the Rydberg states 
results only in single-particle losses [S22], we expect that 
glassy phases should be observable on a time scale suffi¬ 
ciently larger than t and l/7eff. 

These aspects should considerably simplify the detec¬ 
tion of glassy phases in experiments. 


Quench in V/t at fixed T 


Since in experiments it should be easier to perform 
a quench in V/t (which, e.g., with Rydberg atoms 
would entail modifying the laser parameters for Rydberg 
dressing, changing P, or for the depth of the confinement 
to the optical lattice, changing t) than in T, we have 
investigated whether the glassy behaviour described 
above is also found upon a quench in the interaction 
strength. 


We find that similar behavior of and siS^/N to 
the one above is obtained when, starting from a value at 
which the system is superfluid, V/t is abruptly increased 
to an intermediate or large value (see Fig. S4). Following 
this simulation protocol at fixed (low enough) T/t (i.e., 
a quench in V/t) we find no crystalline order and inho¬ 
mogeneity in the thermodynamic limit. Specifically, the 
system is a superglass or a normal glass for intermediate 
or large values oiV/t, respectively. This is entirely anal¬ 
ogous to what found in the main text and above for a 
quench in T. 
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FIG. S4: (color online). Same observables as in Fig. S2 when 
the simulation protocol is based on a quench in the interaction 
strength (see text). Here p = 13/36, N = 2304 and T = 
1/12. The hnal valne of the interaction strength is V/t = 5.4 
[panels (a) and (b)] and V/t = 10 [panels (c) and (d)]. In 
the thermodynamic limit the former choice oiV/t leads to a 
superglass, the latter, to a normal glass. 
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